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Abstract We present independent observations of the solar-cycle variation of flows 
near the solar surface and at a depth of about 60 Mm, in the latitude range ±45°. 
We show that the time-varying components of the meridional flow at these two 
depths have opposite sign, while the time- varying components of the zonal flow are 
in phase. This is in agreement with previous results. We then investigate whether the 
observations are consistent with a theoretical model of solar-cycle dependent merid- 
ional circulation based on a flux-transport dynamo combined with a geostrophic flow 
caused by increased radiative loss in the active region belt (the only existing quan- 
titative model). We find that the model and the data are in qualitative agreement, 
although the amplitude of the solar-cycle variation of the meridional flow at 60 Mm 
is underestimated by the model. 
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1. Introduction 

Solar oscillations are a unique tool to infer conditions inside the Sun. They have been 
recorded with extreme precision since 1996 with the Michelson Doppler Imager (MDI, 
Scherrer et al., 1995) onboard the Solar and Heliospheric (SOHO) space telescope. 
Large-scale rotation inside the Sun can be estimated by inversion of the frequencies 
of millions of global modes of oscillation (e.g., Schou et al., 1998). Rotation is known 
to vary with time in the solar interior at the level of about ±10 ms _1 (e.ij., Schou, 
1999; Howe et al, 2000; Vorontsov et al, 2002; Howe et al, 2006a; Howe et al, 
2006b). These variations, known as torsional oscillations, consist of bands of faster 
and slower rotation that migrate in latitude as the eleven-year solar magnetic cycle 
develops. Torsional oscillations may be driven by the Lorentz force due to a dynamo 
wave (Schiissler, 1981; Yoshimura, 1981; Covas et al, 2000). Other explanations have 
been proposed (see Shibahashi, 2004 and references therein), including the suggestion 
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by Spruit (2003) that torsional oscillations are driven by horizontal pressure gradients 
caused by photospheric magnetic activity. 

Long time averages of surface Doppler measurements have shown the existence 
of a flow from the equator to the poles with an amplitude of 10-20 ms . An 
introduction to the theory of solar meridional circulation is provided by Shibahashi 
(2007). Temporal variations in the meridional flow have been reported by several 
authors. By tracking the small photospheric magnetic features, Komm, Howard, and 
Harvey (1993b), Komm (1994), and Meunier (1999) found a significant change in the 
meridional flow near sunspot latitudes, implying a solar-cycle variation. Variations 
in the surface Doppler meridional velocity have been detected by Ulrich and Boyden 
(2005), in particular at latitudes above 60°. 

Local helioseismology (see, e.g., Gizon and Birch, 2005) also provides reasonable 
measurements of the meridional circulation for latitudes below about 50° (Giles et 
al., 1997). The time-varying component of the meridional flow with respect to a 
long-term average does not exceed ±5 ms" and is consistent with a small near- 
surface inflow toward active latitudes (Basu and Antia, 2003; Gizon, 2004; Zhao and 
Kosovichev, 2004; Gonzalez Hernandez et al., 2006; Komm, Howe, and Hill, 2006) 
and an outflow from active latitudes at depths greater than 20 Mm (Chou and Dai, 
2001; Beck, Gizon, and Duvall, 2002; Chou and Ladenkov, 2001). As discussed by 
Gizon (2004), these variations would appear to be caused by localized flows around 
localized regions of magnetic activity (Gizon, Duvall, and Larsen, 2001; Haber et al., 
2001; Haber et al, 2004). 

The purpose of this paper is two-fold. First, we provide independent measurements 
of the temporal variations of the meridional circulation near the solar surface and 
at a depth of about 60 Mm. For the near-surface meridional flow, we use an original 
technique, which consists of measuring the advection of the supergranulation pattern 
(Gizon, Duvall, and Schou, 2003). The deeper meridional flow is calibrated from 
an earlier time-distance helioseismology observation by Beck, Gizon, and Duvall 
(2002). The meridional circulation measurements at these two depths are compared, 
in particular by looking at the eleven-year periodicity of the flows (Section 2). Our 
study of meridional circulation confirms previous observations (listed above) and is 
complementary to the analysis of zonal flows by Howe et al. (2006a). 

Second, we wish to compare the observations with a theoretical model (Rem- 
pel, 2005; Rempel, Dikpati, and MacGregor, 2005; Rempel, 2006) based on a flux- 
transport dynamo combined with a geostrophic flow caused by increased radiative 
loss in the active region belt, according to Spruit's (2003) original idea. Section 3 
provides a description of the model with a focus on meridional flows, since this aspect 
of the model has not been discussed elsewhere in detail. To our knowledge, this model 
is the only existing quantitative model of the solar-cycle dependence of internal flows: 
it is natural to ask whether it is consistent with the observations. The comparison 
between the observations and the model (Section 4) is encouraging, although some 
inconsistencies cannot be ignored. 

2. Observations of the Meridional Flow 

2.1. Near-Surface Layers 

The method we employ to infer flows near the solar surface is based on the analysis 
of Gizon, Duvall, and Schou (2003), which was originally applied to a single data 
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Figure 1. (a) Rotational velocity (v x ) and (b) meridional flow (v y ) near the solar surface as 
a function of latitude (A). Each MDI dynamics run is plotted with a different color from blue 
in 1996 to red in 2002. The rotational velocity is given with respect to the rotational velocity 
of the small magnetic features (Komm, Howard, and Harvey, 1993a). 



set from 1996. Here we use a series of MDI full-disk Doppler images covering the 
period 1996-2002. Each year, two to three months of continuous Dopplergrams are 
available for analysis (MDI Dynamics runs). The MDI data after 2002 were not used 
simply because we are analyzing an existing pre-processed data set. 

Dopplergrams were tracked at the Carrington rotational velocity to remove the 
main component of rotation. We used /-mode time-distance helioseismology (Duvall 
and Gizon, 2000) to obtain every 12 hour a 120° x 120° map of the horizontal 
divergence of the flow field 1 Mm below the photosphere. The main component of 
the divergence signal is due to supergranulation. For any given target latitude (A) 
we considered a longitudinal section of the data 10° wide in latitude. Using a local 
plane-parallel approximation in the neighborhood of latitude A, the divergence signal 
was interpolated onto a Cartesian grid sampled at 2.92 Mm in the x (prograde) and 
y (northward) coordinates. The divergence signal was decomposed into its harmonic 
components exp(i(fc x a; + k y y — Lot)) to obtain a local power spectrum P(k,u>; A), 
where k = (k x ,k y ) is the horizontal wavevector and ui is the angular frequency. At 
fixed fci?0 = 120, we fit for two functions / and g and a horizontal vector v such that 
P(k, lo; A) = f{k) g(u — k ■ v(A)). This representation fits the data adequately. As 
was done by Gizon, Duvall, and Schou (2003), we interpret v(X) to be a horizontal 
flow causing a Doppler shift Alu = k ■ v. This flow is likely to be an average over 
the supergranulation layer, which has been estimated to reach depths greater than 
10 Mm by Zhao and Kosovichev (2003). 

Figure 1 is a plot of v x (X) and v y (X) for each full-disk MDI run as a function of 
latitude in the range | A | < 50° . To reduce random noise, the North - South symmetric 
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component of v x and the antisymmetric component of v y are extracted. Over the 
period 1996-2002, v x varies by 12 ms _1 peak-to-peak at the Equator (Figure la). 
The Meridional flow is poleward with a mean amplitude of 10 ms _ at latitude 20° 
(Figure lb). The peak-to-peak variation of the meridional flow is 7 ms at A = 30°, 
i.e. a significant fraction of the time-average value. We estimate that the standard 
deviation of the noise at a particular latitude (5° bin) for any given year is less than 
1 ms~ . The systematic errors that depend on position on the solar disk have been 
measured to be very low (less than 5 ms _1 over the 120° x 120° region of analysis). 

2.2. Deeper Inside the Sun 

In order to probe deeper layers into the solar convection zone, we used acoustic waves 
and time-distance helioseismology. For each three-month period, travel times were 
measured by cross-correlation of the Doppler oscillation signal recorded during the 
MDI structure program (nearly continuous coverage but lower spatial resolution) 
according to the procedure described by Giles (1999). Using a mean travel distance 
of 17° enables us to probe layers about 60 Mm below the surface. The full details of 
this analysis can be found in Beck, Gizon, and Duvall (2002). Waves that propagate 
in the North - South direction are used to infer the meridional flow, while waves that 
propagate East -West are used to infer the zonal flows. In order to convert travel- 
time shifts into flows in units of ms , we use a simple calibration of v x at a depth 
based on the observation by Howe et al. (2006a) (global-mode helioseismology) that 
the amplitude of the time- varying component of the zonal flow is nearly independent 
of depth. We choose the near-surface zonal-flow measurements of Section 2.1 as 
a reference. The calibration of v x is then used to calibrate v y . We find that the 
meridional flow at a depth of 60 Mm is poleward at all latitudes and has a maximum 
value of 6 ms _1 at latitude 25°. For a particular year and at fixed latitude (5° bin), 
the standard deviation of the noise is about 2 ms -1 , significantly more than for the 
near-surface measurements. 

2.3. Solar-Cycle Variations 

In order to discuss the solar-cycle dependence of the flows and to study the phase 
relationship between the flows measured at the two different depths, we extract the 
eleven-year periodic component from the data, as was done by Vorontsov et al. (2002) 
and Howe et al. (2006a) for zonal flows. At each latitude A and for each depth, we 
fit a function of the form 



to the observed velocity Vi(X,t), where the index i refers to either the x or the y 
component of the flow. The long-term average is given by Vi, while the amplitude and 
the phase of the time-varying component are denoted by v[ and <f>i respectively. We 
extract a 11-year periodicity from the data, since it is known from other observations 
that this is the dominant mode (we do not determine the 11-year periodicity based 
on the dataset itself). Shorter and longer periodicities are certainly present in the 
data; however, the length of the dataset does not allow for a determination of the 
full spectrum of modes. 




(1) 
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Figure 2. Eleven-year periodic component of the meridional (top two panels) and zonal 
(bottom two panels) flows as function of time and latitude at two different depths in the 
solar interior: near the surface (top and third panels) and 60 Mm deep (second and bottom 
panel). The color bar is in units of ms — . A positive value indicates a poleward meridional 
or prograde zonal flow. The observations (ti; — Vi) cover the first six years, while the purely 
sinusoidal component (vi — Vi) is extrapolated in time (beyond the white vertical white line). 
The black curves indicate the mean latitude of magnetic activity. 
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Figure 3. Amplitude (i/) of the eleven-year periodic component of the (a) meridional and 
(b) zonal flows. The near-surface values (thick solid lines) are absolute measurements. The 
calibration of the observations at 60 Mm depth (thin lines) follows from the assumption that 
the amplitude of the zonal torsional oscillation (panel b) is independent of depth over the 
latitude range |A| < 45°. 




Figure 4. Phase difference [A</> = </>(deep) — </>(surface)] between the eleven-year periodic 
components of the flows measured at a depth of 60 Mm and near the surface. The solid line is 
for the meridional flow and the dashed line is for the zonal flow. 
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The eleven-year periodic components of the meridional and zonal flows are shown 
in Figure 2. The torsional-oscillation pattern is clearly seen at both depths with 
an amplitude and a phase comparable to previous measurements (e.g., Howe et ai, 
2006a). The meridional flow also contains a significant eleven- year periodic compo- 
nent. Near the solar surface, the residuals indicate the presence of a North -South 
inflow toward the mean latitude of activity (e.g., Zhao and Kosovichev, 2004; Komm, 
Howe, and Hill, 2006), while the data are consistent with a horizontal outflow from 
the mean latitude of activity deeper into the convection zone (e.g., Chou and Dai, 
2001; Beck, Gizon, and Duvall, 2002). 

Figure 3 gives the amplitudes of the time- varying components of the flows (v'i). 
Under the assumption (Section 2.2) that v x does not vary appreciably with depth 
(5 ms -1 latitudinal average), then the amplitude of the time- varying meridional flow 
(v' y ) is also found to be approximately independent of depth (v' y ~ 3 ms — at 20° 
latitude) . The evidence that the time- varying components of the meridional flow near 
the surface and deeper in the interior are anti-correlated is given in Figure 4, which 
shows the difference in <j) y at the two depths. On the contrary, there is no significant 
phase variation with depth for the zonal flow. 

3. Theoretical Model of Time- Varying Flows 

The model results presented here are based on a non-kinematic flux-transport dy- 
namo model developed recently by Rempel. This model combines the differential 
rotation and meridional flow model of Rempel (2005) with a flux-transport dynamo 
similar to the models of Dikpati and Charbonneau (1999) and Dikpati and Gilman 
(2001). We emphasize that this model is intended to give a fundamental understand- 
ing of the basic cycle properties and their relation to observable variations of zonal 
and meridional flows. Therefore we focus here only on axisymmetric and North - 
South averaged quantities. Details of the model can be found in Rempel (2006). Since 
a detailed comparison with observed torsional oscillations can be found in Rempel 
(2006) and Howe et al. (2006b), we focus here on the meridional flow variations. 

The differential rotation model utilizes a meanfield Reynolds-stress approach that 
parametrizes the turbulent angular momentum transport (Kitchatinov, Riidiger, 
1993; A-effect) leading to the observed equatorial acceleration. In this model the 
tachocline is forced through a uniform rotation boundary condition at the lower 
boundary of the computational domain. A meridional circulation, as required for a 
flux-transport dynamo, follows self-consistently through the Coriolis force resulting 
from the differential rotation. 

The computed differential rotation and meridional flow are used to advance the 
magnetic field in the flux-transport dynamo model, while the magnetic field is allowed 
to feed back through the meanfield Lorentz-force (J) x (B) (the contribution of the 
fluctuating part (J' x B') is not well known and neglected here). 

We find in our model that the Lorentz-force feedback can only account for the 
poleward propagating branch of the torsional oscillations, while the equatorward 
propagating branch in latitudes beneath 30° requires additional physics. Parametriz- 
ing the idea proposed by Spruit (2003) that the low-latitude torsional oscillation is a 
geostrophic flow caused by increased radiative loss in the active region belt (due to 
small scale magnetic flux) leads in our model to a surface oscillations pattern in good 
agreement with observations. In order to force a torsional oscillation with around 1 
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nHz amplitude, a temperature variation of around 0.2 K is required. As a side effect 
the cooling produces close to the surface (in our model at r = 0.985 Rq) an inflow 
into the active-region belt of around 2.3 ms . 

We incorporated this process by adding a surface-cooling term that is dependent 
on the toroidal-field strength at the base of the convection zone, which is assumed 
to be the source for active-region magnetic field (the small scale flux required for the 
surface-cooling is a consequence of the decay of active-regions). Observations show 
that the low-latitude branch of torsional-oscillations starts around one to two years 
before the sunspots of the new cycle appear, ft is possible that magnetic flux rises 
towards the surface without forming sunspots in the beginning of a cycle providing 
enough small-scale magnetic field, this is however currently neither confirmed nor 
ruled out by observations. Alternative explanations for the low-latitude branch of 
torsional oscillations such as the models of Sehussler (1981), Yoshimura (1981), and 
Covas et al. (2000) are based on the longitudinal component of the Lorentz force. 
Recently (Rempel, 2007) showed that torsional-oscillations forced that way are close 
to the Taylor-Proudman state (alignment of phase with the axis of rotation), which 
contradicts observations. In addition, the resulting meridional surface-flow has the 
wrong sign (active region belt outflow). Despite some shortcomings, the model of 
Spruit (2003) is currently the only proposed explanation that is consistent with the 
observed meridional and zonal-flow variations close to the solar surface. 

Figure 5 summarizes the results of the model in latitudes below 45°. Figure 5a 
shows the temperature fluctuation (color shades) caused by increased surface cooling 
in the active region belt. The contour lines indicate the magnetic butterfly diagram 
computed from the toroidal field at the base of the convection zone in the model. 
Figure 5b shows cycle variations of the angular velocity (torsional oscillations) and 
Figure 5c shows cycle variations of the horizontal meridional flow velocity. At the 
equatorward side of the active region belt (indicated by the butterfly diagram) the 
rotation rate is increased, which is consistent with the increased poleward meridional 
flow transporting material toward the axis of rotation. On the poleward side of the 
active region belt the rotation rate is lower, while the meridional-flow perturbation 
is equatorward. At a depth of around 50 Mm (Figure 5d) the meridional-flow per- 
turbation is almost anti-correlated to the surface flow (active-region belt outflow), 
indicating that the surface cooling drives a flow system that closes in the upper third 
of the convection zone. The flow amplitude at a depth of 50 Mm is around one order 
of magnitude lower compared to the surface flow due to the significant increase in 
density. 

4. Discussion 

We presented a model that is qualitatively consistent with the observations, in partic- 
ular the phase of the solar-cycle variations of the flows. Since the surface-cooling term 
is parametrized, the model can only make predictions about the relative amplitude 
of zonal and meridional flows at different depths, and not about the absolute values. 
Near the surface, the model is in agreement with the data: the relative amplitudes 
of the torsional oscillation and the time-varying component of the meridional flow 
and their relative phase are reproduced well. Deeper in the interior, it appears that 
the model underestimates the amplitude of the time variations of the meridional 
flow by an order of magnitude; however, the flow variation is in antiphase to the 
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Figure 5. Model results, (a) Surface temperature variation (blue: cold, red: hot, amplitude: 
0.2 K). (b) Torsional oscillations (blue: slower, red: faster rotation, amplitude: 1.35 nHz). (c) 
Meridional flow variation at r = 0.985 iJ© (blue: equatorward, red: poleward motion, amplitude: 
2.3 ms -1 . (d) Meridional flow variation at r = 0.93 Rq (blue: equatorward, red: poleward 
motion, amplitude: 0.22 ms -1 ). The variation of the meridional flow pattern at r = 0.985 Rq 
is almost in anti-correlation with the flow at r = 0.93 Rq (~ 50 Mm depth). In all four panels 
the contour lines indicate the butterfly diagram computed from the toroidal field at the base 
of the convection zone. 
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surface flow as seen in the data. The lower velocity at depth in the dynamo model 
is a consequence of mass conservation (strong increase in density with depth). The 
much larger outflow which is observed in the data cannot be balanced by an inflow 
close to the surface unless the outflow is confined to a very narrow layer. 

Overall, it is fair to say that the model is encouraging. Alternative models ex- 
plaining low-latitude torsional oscillations through the action of the longitudinal 
Lorentz force tend to produce zonal flow patterns in the Taylor-Proudman state 
(Rempel, 2007) and meridional outflows rather than inflows close to the surface. On 
the other hand, zonal flows appear during solar minimum when no active regions 
are present. This requires, in the case of thermal forcing, that a sufficient amount 
of magnetic flux is present in the form of small scale flux elements not evident in 
synoptic magnetograms (see e.g., Section 6.3 Spruit, 2003). Also a local treatment 
of the regions of strong magnetic-field concentrations (sunspots and active regions) 
might be necessary to obtain a better match between the model and the data. 

On the observational side, it would be useful to invert the travel-time measure- 
ments in order to obtain improved and more reliable estimates of the depth variations 
of the flows. 
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